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Abstract. 

The GALPROP cosmic -ray (CR) propagation model has 
been extended to three dimensions including the effects of 
stochastic SNR sources, a comprehensive cross-section data- 
base, and nuclear reaction networks. A brief description of 
the new code is given and some illustrative results for the 
distribution of CR protons presented. Results for electrons 
and 7-rays are given in an accompanying paper. 



1 Introduction 

We have previously described a numerical model for the 
Galaxy encompassing primary and secondary cosmic rays, 
7-rays and synchrotron radiation in a common framework 
(Strong et al. 2000 and references therein). Up to recently 
our GALPROP code handled 2 spatial dimensions, {R, z), 
together with particle momentum p. This was used as the 
basis for studies of CR reacceleration, the size of the halo, 
positrons, antiprotons, dark matter and the interpretation of 
diffuse continuum 7-rays. 

Some aspects of the problem cannot be addressed in such 
a cylindrically symmetric model: for example the stochastic 
nature of the cosmic-ray sources in space and time, which is 
important for high-energy electrons with short cooling times, 
and local inhomogeneities in the gas density which can affect 
radioactive secondary /primary ratios. 

In common with most other models it has been previously 
assumed that the CR source function can be taken as smooth 
and time-independent, an approximation justified by the long 
residence time (> 10^ years) of cosmic -rays in the Galaxy. 
However the inhomogeneities have observable consequences, 
and their inclusion is a step towards of the goal of a "re- 
alistic" propagation model based on Galactic structure and 
plausible source properties. The original motivation for this 
extension was to study the high-energy electrons, since the 
observation of the > 1 GeV excess in the EGRET spec- 
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trum of the Galactic emission has been proposed to origi- 
nate in inverse-Compton emission from a hard electron spec- 
trum; this hypothesis can only be reconciled with the local 
directly-observed steep electron spectrum if there are large 
spatial variations which make the spectrum in our local re- 
gion unrepresentative of the large-scale average. 

Here we briefly describe an extension of the model to 3D, 
which can address these issues, and illustrate the results for 
protons. The new network and cross-sections is illustrated 
for B/C in the 2D case. The effect on electrons and 7-rays is 
presented in an accompanying paper (Strong and Moskalenko, 
'A 3D time-dependent model for cosmic rays and 7-rays', 
these proceedings: hereinafter paper II). 

2 Model 

The GALPROP code, which solves the CR propagation equa- 
tions on a grid, has been entirely rewritten in C++ using the 
experience gained from the original version and including 
both 2D and 3D spatial grid options. The 2D mode essen- 
tially duplicates the original version, with improved cross- 
section routines. In 3D [x, y, z,p) the propagation is solved 
as before using a Crank-Nicolson scheme. The additional 
dimension considerably increases the computer resources re- 
quired, but a 200 pc grid cell or finer is still practicable. As 
in the original version, the effects of diffusion, convection, 
diffusive reacceleration, and energy losses are included, each 
with adjustable parameters defined for a GALPROP run. 

2. 1 Cross-sections and reaction network 

Cosmic-ray nuclear reaction networks are included with a 
comprehensive new cross-section database; this allows the 
models to be tuned on stable and radioactive CR second- 
ary/primary ratios, in particular B/C and ^"Be/^Be. The nu- 
clear reaction network is built using the Nuclear Data She- 
ets. The isotopic cross section database consists of more 
than 2000 points collected from sources published in 1969- 
1999. This includes a critical re-evaluation of some data and 
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cross checks. The isotopic cross sections for B/C were calcu- 
lated using the authors' fits to major beryllium and boron pro- 
duction cross sections C,N,0 Be,B- Other cross sections 
are calculated using the semi-phenomenological approxima- 
tions by Webber et al. (1990) (code WNEWTR.FOR of 1993) 
and/or Silberberg and Tsao (code of 2000) renormalized to 
the data where it exists. 

The reaction network is solved starting at the heaviest nu- 
clei (i.e. ^^Ni). The propagation equation is solved, com- 
puting all the resulting secondary source functions, and then 
proceeds to the nuclei with A — 1. The procedure is repeated 
down to A = 1. In this way all secondary, tertiary etc. reac- 
tions are automatically accounted for. To be completely ac- 
curate for all isotopes, e.g. for some rare cases of /3^-decay, 
the whole loop is repeated twice. Our preliminary results 
for all cosmic ray species < 28 are given in Strong and 
Moskalenko (2001). 



2.2 Stochastic sources 



Another major enhancement is the inclusion of stochastic 
SNR events as sources of cosmic rays. The SNR are char- 
acterized by the mean time Isnr between events in a 1 kpc^ 
unit volume, and the time tcR during which an SNR actively 
produces CR. 

The propagation is first carried out for a smooth distri- 
bution of sources to obtain the long timescale solution, us- 
ing the fast technique described in Strong and Moskalenko 
(1998); then the stochastic sources are started and propa- 
gation followed for the last 10^ years or so with timesteps 
sufficiently fine (e.g., 10"^ years) to resolve the SNR events 
and all propagation effects (see below for more details of the 
method). For high-energy (TeV) electrons which lose energy 
on timescales of 10^ years the effect is a very inhomoge- 
neous distribution with consequences for diffuse 7-rays, as 
shown in the accompanying paper 11. However also the pro- 
tons (and other nuclei) show fluctuations which are strongly 
energy-dependent. The amplitude of the fluctuations of both 
electrons and nuclei depends on the parameters tsNR and 
tcR- The time tsNR is adjusted to be consistent with esti- 
mates of the SNR rate (e.g., Dragicevich et al. 1999); models 
for shock acceleration in SNR indicate 10^ < tcR < 10^ yr, 
the sources switching off when they move from the adiabatic 
to the radiative phase (Stumer et al., 1997). 

An additional advantage of the 3 spatial dimensions is that 
the gas distribution (for energy losses, secondary production 
and 7-rays) can be modelled in as much detail as required, 
based on current HI, CO, FIR, etc. surveys. Spiral struc- 
ture and local inhomogeneities can be included if required. 
Global parameters such as the CR luminosity of the Galaxy, 
and hence the average SNR energy injection into CR, are also 
computed since they provide essential constraints (see paper 

n). 



2.3 Extension to 3D 

Since the solution in {x, y,z,p) on a fine grid involves large 
arrays and many time-steps, the code has been made vec- 
torizable so that it can benefit from the use of vector ma- 
chines with a speed gain typicaUy 100. This is only required 
for the final "stochastic SNR" part of the solution, since the 
"smooth" solution is fast enough on non-vector machines. 
In addition we have the option to make use of symmetry in 
the spatial dimension where it does not affect the application 
of the results (e.g., about z = 0). This leads to an addi- 
tion large gain in speed. Typical parameters are: grid size 
Ax = Ay = Az = 200 pc. Galaxy dimension 60 kpc x 60 
kpc X 4 kpc, and energy range 1 MeV - 1 TeV on a logarith- 
mic scale with factor 1.2. ^ 



3 CR Protons 

For illustration we show results for a model with reaccelera- 
tion based on Strong et al. (2000), and tsNR = 10^ yr corre- 
sponding to a Galactic rate of 3 SN/century. Fig. 1 shows the 
distribution of protons in the Galactic plane (z = 0) for a rep- 
resentative quadrant, at various energies. The stochastic SNR 
source produce fluctuations, which are a minimum around 
1 GeV and increase at low energies due to energy losses and 
at high energies where the storage of particles in the Galaxy 
is much reduced so that the effect of sources manifests itself 
on the distribution. Note that the nature of the fluctuations is 
different at low and high energies. 

Fig. 2 shows a sample of proton spectra at different {x, y) 
positions; the fluctuations are evident but much smaller than 
for electrons (paper II). 

The fluctuations in the GeV nuclei will have some effect 
on the 7r° -decay diffuse 7-ray emission above 100 MeV, but 
the variations will be small compared to the inverse Compton 
component since the GeV nuclei variations are much smaller 
than for the TeV electrons responsible for GeV 7-rays via IC 
(see paper II). Above 100 GeV 7-ray energies however the 
effect on the 7r°-decay component wiU be larger, and should 
be detectable for example by the GLAST mission. 



4 Secondary/primary ratio 

We illustrate the new cross-sections and network with the 
B/C ratio (Fig. 3) for the case of the same model (with reac- 
celeration) as for electrons. More details of the application 
to nuclei can be found in Strong and Moskalenko (2001) 
and Moskalenko et al. 'New calculation of radioactive secon- 
daries in cosmic rays' (these proceedings) and to antiprotons 
in Moskalenko et al. 'Secondary antiprotons in cosmic rays' 
(these proceedings). 

' As usual the software will be made available on the WWW at 
http://www.gamma.mpe-garching.mpg.de/~aws/aws.html 




Fig. 1. Flux of protons with E= 64 MeV, 1 GeV, 8.2 GeV, 130 GeV, and 1 TeV (left to right, top to bottom) at 2 = 0, for tsNR = 10* yr. 
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Fig. 2. Sample proton spectra for various x,y at z — 0, illustrating 
the variations of shape with position. Model as Fig. 1. 



Fig. 3. B/C ratio as calculated for a model with reacceleration using 
the new reaction network and cross-sections. Lower curve: inter- 
stellar, upper curve: modulated for 450 MV. Data: see Strong and 
Moskalenko (2001). 



5 Conclusion 

This paper is intended only as an illustration of the possi- 
bilities opened up with such a 3D code. Future work will 
concentrate on the effect of inhomogeneities in the gas dis- 
tribution (e.g., the local bubble) on radioactive CR nuclei, as 
well as impUcations for electrons and 7-rays. 
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